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Abstract. Temperature dependent photoemission spectroscopy in Lio.gMogOiy 
contributes to evidence for one dimensional physics that is unusually robust. Three 
generic characteristics of the Luttinger liquid are observed, power law behavior of the 
k-integrated spectral function down to temperatures just above the superconducting 
transition, k-resolved lineshapes that show holon and spinon features, and quantum 
critical (QC) scaling in the lineshapes. Departures of the lineshapes and the scaling 
from expectations in the Tomonaga Luttinger model can be partially described by 
a phenomenological momentum broadening that is presented and discussed. The 
possibility that some form of Id physics obtains even down to the superconducting 
transition temperature is assessed. 
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1. Introduction 

Lio.gMo60i7, the so-called lithium purple bronze (hereafter called LiPB), has emerged as 
a paradigm quasi one-dimensional (Id) material for studying the stability of Luttinger 
liquid chains against single particle electron hopping t± between the chains p] . The usual 
theoretical understanding [2j [3] and empirical experience is that quasi-ld materials can 
display such Id behavior at sufficiently elevated temperatures but that, with decreasing 
temperature (T), t± leads to a crossover to some kind of 3d behavior, e.g., a Fermi 
liquid, or an ordered state such as a charge density wave (CDW) or spin density wave 
(SDW). Because of the excellent Fermi surface (FS) nesting that is typical of quasi-ld 
materials, density waves are especially likely. Focusing on photoemission spectroscopy 
but considering other data as well, this paper explores the possibility that LiPB, in spite 
of excellent FS nesting, does not develop a density wave and crosses over to 3d behavior 
only when it becomes a superconductor (SC). If such is true then the SC could well be 
unconventional, as suggested in a recent paper jl] reporting that the superconducting 
critical field is far above the Pauli limit. 

LiPB is a material with overall 3d bonding but whose valence band electronic 
structure is derived from two parallel zig-zag Mo-0 chains per unit cell. The chains 
lie in well separated planes and receive electrons donated by Li ions located out of the 
planes. Early tight binding [5] and more recent LDA [6] band calculations agree that the 
Mo-0 orbitals of the chains give rise to four bands, two of which (A, B) are always below 
the Fermi energy (E-p) and two of which (C,D) cross Ep. The resulting FS is quasi-ld 
but the LDA calculations also show some splitting and warping due to single particle 
hopping £j_ between chains both within a unit cell and in neighboring unit cells. Angle 
resolved photoemission spectroscopy (ARPES) observes all four bands with the general 
behavior that is predicted in LDA calculations [7J El [9], including the quasi-ld FS [7J. 
Figures presented in later discussions of ARPES data show the four bands (Fig. E]Ja) in 
Section 3), and the FS and Brillouin zone (Fig. @|c) in Section 4). The magnitude of t± 
is discussed in Section 6 of the paper. 

The physical properties of LiPB were first investigated in the 1980 's [101 [TTJ [121 H3J 
HU [15j [16j [TTJ, HU [T9l [20] . Electrical transport [10] was found to be highly anisotropic and 
superconductivity (SC) [TOj [HI [151 HB1 HE] was observed below Tsc = 1-9 K, albeit not 
in all samples [19] and with a variation of Tsc [El El IISI [El [19] that might involve the 
Li stoichiometry [21] or disorder [15]. More recent studies [221 E31 EH [251 EEl [27] confirm 
the anisotropic transport but with differences in the absolute values of the resistivities 
along various axes, and also in the magnitudes of the anisotropy ratios. These differences 
may reflect variation in samples and in the method [23] of the measurement. 

2. Resistivity Upturn and Evidence Against a Density Wave 

The early experimental work found that the resistivity is metallic with decreasing T 
only down to a minimum at T m ; n fa 26 K, below which there is an abrupt upturn. This 
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Figure 1. Panel (a) shows the temperature-dependent resistivity as taken from the 
publication of Xu et al. [33] together with a line-fit of the data with the function 
p Flt = ATP + BT~>. The coefficients found were (3 = 0.53 ± 0.1 and 7 = -1.54 ± 0.2. 
The inset shows that the data from Ref. [TU] and Ref. [25] are qualitatively similar. 
Note that the inset shows the resistivity normalized to its 300K value in order to be 
comparable with that of Ref. [25J. Panel (b) shows the gap as extracted in Ref. [21)] 
together with the gap extracted here by using our fit Apu = 2 fcsT ln(pFit/ Po)- Also 
included is a mean-field temperature behavior of a gap A m f — 3.52 fe^Ic yl — T/Tc 
calculated with Tc=26K. Note that for the much larger mean-field gap, the axis is on 
the right. 

finding has been robust over time. Fig. Q^a) shows early data [TD] and more recent data 
[25] [26] for the upturn as measured along the Id direction. An upturn at the same 
temperature is also seen for the two directions perpendicular to the chains. There is 
yet no consensus as to the origin of this upturn because each proposed explanation 
requires that one or another piece of data be ignored, be given less weight in the 
argument, or be left unexplained [28]. Early discussions ascribed the upturn to Anderson 
localization due to disorder [121 1131 EH ES] or to a single particle gap, most likely due 
to CDW formation [lOl ETJ EH EHJ [20] . Both direct and nuanced lines of argument 
[IQI HU [121 [TU [131 HH [191 En! [22I [23I [2II [25l EHl I2Z] for and against these two hypotheses 
have been presented many times in the transport literature. Ref. [22] is the most 
comprehensive recent overview of the issue. 

The disorder hypothesis was proposed [121 El] because a number of important 
expected signatures of a density wave are missing in LiPB. X-ray diffraction studies 
that have revealed CDW ordering in related materials [29] have not observed a CDW 
in LiPB in spite of repeated attempts over time |30j. Low-T optical spectroscopy that 
has observed single particle gaps in various CDW and SDW materials [TF1 EH [32] finds 
[T71 [22] in LiPB no gap even down to 6 K and 1 meV. This finding would seem to be 
conclusive for the magnitude of gap (~ 8 meV) expected in a mean field description of 
a CDW transition at T min , as seen from the BCS curve of Fig. E^b). Comparison to 
mean field is done here for lack of a better alternative. On the one hand, the actual 
transition temperature for a quasi- Id system is likely controlled by a small t± and so 
can be much less than the mean field temperature for a single chain [33]. In this case 
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the actual gap could be much larger than that obtained in this estimate. On the other 
hand the single chain gap for T=0 can be much reduced from its mean field value by LL 
fluctuations on the chain, in which case the gap could be smaller than in this estimate. 
In fact Fig. []Jb) shows that the T-dependent gap that is deduced [26] by describing the 
upturn as an Arrhenius law is slightly smaller than 1 meV at low T and so could have 
gone undetected in the optical work. Such a small gap was invoked as an explanation for 
the recent and very exciting finding [26] that a sufficiently large magnetic field applied 
along the Id direction can entirely suppress the upturn. The explanation put forth is 
that the magnetic field quenches the CDW and hence the gap because Zeeman splitting 
of the bands disrupts FS nesting. The discussion in Ref. [26] also took up a previous 
proposal [21] for a purely electronic CDW, driven entirely by Coulomb interactions in 
a strong coupling limit. This proposal was made to overcome indications from optical 
spectroscopy [22] and thermal expansion [21] measurements that the lattice has no 
strong involvement with the upturn, which implies that a conventional CDW involving 
a periodic lattice distortion does not occur. While such a purely electronic CDW would 
be more difficult to detect in the conventional x-ray diffraction performed for LiPB to 
date [30] , it is certainly not impossible as shown by such an observation for a Bechgard 
salt [31]. 

The great beauty of the experimental magnetoresistance result [26J and the neatness 
of the explanation notwithstanding, there is no CDW theory that would explain such 
an unusual gap function and one must also ignore another important missing signature 
of a single particle gap [lOj [T21 E2], that the dc magnetic susceptibility measured by 
three groups over time is entirely unaffected by the resistivity upturn at T min . Ref. 
[22] points out that because the magnitude of the measured susceptibility involves a 
near cancelation of several large negative and positive contributions the measurement 
is particularly sensitive to a change in any one of them, e.g that of the conduction 
electrons. This is perhaps the strongest single piece of all the many evidences that 
collectively weigh heavily against a CDW. The lack of CDW signatures in LiPB is all 
the more striking because the closely related compounds KM06O17 and NaMoeOi7 do 
exhibit CDW behavior clearly seen with x-ray diffraction [29] and optics [17] . These 
materials are also planar but have higher symmetry such that the planes contain three 
Mo-0 chains oriented at 120 degrees to one another. The CDW Q- vector is such as to gap 
out the quasi- Id bands associated with two of the three chains. Finally the possibility 
of an SDW in LiPB has been excluded by muon spin relaxation measurements [35] . The 
last section of the paper returns to the issue of the upturn and to the other parts of 
Fig.E 

3. Signatures of Luttinger liquid behavior 

Early workers did not include the ideas of the Luttinger liquid (LL) in their thinking. LL 
behavior in LiPB was first demonstrated by photoemission experiments. Indeed, these 
experiments were undertaken specifically because of the missing signatures for CDW 



PES and the Unusually Robust Id Physics of Lio gMo^On 



5 



formation. Pioneering studies [35] had observed in angle integrated photoemission of 
certain quasi- Id materials one of the characteristic features expected of LL behavior, 
that the angle integrated single particle spectrum approaches Ep as a power law. These 
studies suffered from some ambiguity because the materials studied all manifested 
low T CDW formation and CDW fluctuations at higher T can in principle cause 
pseudogap behavior [331 [37] that might mimic the expected LL behavior. Thus LiPB 
was conceived [7J as a non-CDW quasi-ld paradigm for studying LL behavior. It was 
also envisioned that ARPES experiments could distinguish features of the single particle 
spectral function that are unique to LL behavior. A series of ARPES studies dating 
from 1999 [Zl El El EE! EH1 SUl SH S3 S3] has born out this motivation. This section 
gives a guide to the key results of these various past ARPES studies and points out 
two other measurements using different techniques that have added very importantly to 
the case for LL behavior. Sections 4 and 5 present new photoemission data and a new 
phenomenological analysis of ARPES data, respectively. 

ARPES measures the single particle spectral function. For the paradigm one 
band Tomonaga-Luttinger (TL) model [44J the spectral function [45] shows three 
characteristic signatures of LL behavior. The ARPES measurements described below 
have demonstrated all three of these signatures. First, as already mentioned, the 
momentum(k)-summed spectrum approaches Ep as a power law characterized by an 
anomalous exponent a. Second, the k- resolved single particle lineshape is composed of 
two features, a holon peak with a leading spinon edge, the two dispersing with different 
velocities, v c , v s , respectively. This lineshape signifies the absence of quasi-particles and 
the fractionalization of the electron into density fluctuation modes of charge (holon) and 
spin (spinon). Third, the TL-model is quantum critical (QC). QC systems display scale 
invariance, i.e. at T=0 K and for large distances and long times, the correlation functions 
lack a characteristic scale and take the simplest possible scale-free, functional form, a 
power law. Departures from T=0 K satisfy simple scaling laws [46l 147] with T the only 
scale. For the case of a spin rotationally invariant interaction the TL-model spectral 
function [45] explicitly obeys the ideal scaling form A(k, u, T) = T v A(vk/T,u/T) where 
A is a universal scaling function, k is measured from the Fermi momentum kp, u is 
measured from the Fermi energy Ep, and v is a constant with units of velocity. Within 
the one band TL-model the T-scaling exponent r\ is then also determined by a scaling 
relation to be r\ = (a — 1). 

The initial ARPES studies [7] were performed at relatively high temperatures of 
200 K to 300 K to avoid any influence of the putative CDW transition at T min but the TL- 
model theoretical lineshapes then available [48] for comparison [7_J [38] to the ARPES 
data were for T=0 K. Over time the quality of the spectra improved and nonzero-T 
theory lineshapes [45] became available. The latter were particularly important because 
including the effect of temperature changed the parameters deduced from the T=0 K 
comparison such as to improve the internal consistency of the description. It was found 
thereby [391 HQ] that the high T ARPES spectra could be generally well described by 
the theory lineshapes for the measurement T and for an a value the same as determined 
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directly from the angle integrated spectrum. These results can be seen in Figs. 7 and 8 
of Ref. [39J and Figs. 4,5 and 6 of Ref. [ID]. Fig. 7 of Ref. [39] sketches the bands and 
shows ARPES data for T=250 K. Fig. 8 of Ref. [39J compares the data to theoretical 
lineshapes calculated for T=250 K, for a=0.9 and for a range of values of ^. The 
lineshapes include the experimental broadening in k and u. The choice of ^=2 gives 
the best agreement with the data. Fig. 6 of Ref. [40J shows a similar comparison of the 
data to theory for a range of values of a and one can see that a=0.9 gives the best 
agreement for the rate of falloff of the ho Ion peak intensity as k approaches kp. Fig. 7 
of Ref. [39] is interesting for showing explicitly that k-integration of ARPES data for a 
Fermi liquid material yields a Fermi edge whereas one obtains a power law at for 
LiPB. 

Fig. 8 of Ref. [39J (same as Fig. 5 of Ref. [40]) makes an important point concerning 
the ARPES lineshapes. For a > 0.5, as is the case for LiPB, the spinon feature of 
the lineshape is an edge singularity rather than a peak singularity. Including also the 
broadening due to temperature and experimental resolutions yields a lineshape which at 
first glance does not have two distinct features and so could arouse skepticism as to the 
claim of observing spin-charge separation. Nonetheless the various panels of the figure 
make it clear that the two features are indeed present and are still very visible because 
the spinon edge disperses at a rate that differs from that of the holon peak, controlled 
by varying ^. Thus these lineshapes are quite unique, are characteristic of LL behavior, 
and are much different from the usual dispersing peaks seen in ARPES. 

Two other important findings for the ARPES lineshapes of LiPB are that the same 
spectra are obtained [41] for samples made using each of the two main crystal growth 
techniques and that, within experimental resolutions, the same spectra are obtained 
when measured [42] at a photon energy hu=500eV as when measured with the lower 
photon energies hp between 20 eV and 30 eV that were used for all the other ARPES 
summarized here. The significance of the former is that it dispels the possibility of 
sample growth method being the origin of an early ARPES report [49j of Fermi liquid 
lineshapes for T above T m in and of an 80 meV gap for T below T m j n . This anomalous 
finding was disputed [501 15T] and since then has never been replicated. The significance 
of the latter is that the higher photon energy spectra are more bulk sensitive and the 
agreement of the spectra is consistent with the likelihood based on the crystal structure 
that the quasi-ld chains lie two layers below the cleavage plane and are hence well 
protected from surface effects, the same situation that has made ARPES relevant to the 
bulk properties of many superconducting cuprate materials. 

Scanning tunneling spectroscopy (STS) [52] has made an important contribution 
to the case for LL behavior. This work is notable for its T-range from 55 K down to 5 
K, only slightly above Tsc- The experimental resolution deduced from the STS spectra 
is 9 meV, essentially the same as the gap value implied by a mean field transition at 
T min . STS also observes a power law density of states at E F , although the values of 
a ~ 0.6 are smaller than those measured in the high temperature ARPES. No change 
was detected in the T-range of the resistivity upturn, consistent with the lack of a gap 
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Figure 2. ARPES spectra of Lio.gMo60i7. (a) Intensity map of 300K spectra of 
T-dependent data set analyzed quantitatively in Ref. |43) . D is shown as a dashed line 
because it is strong HO E] only for different experimental geometries and k-paths. 
(b) Overplot of spectra in range of box in (a). Different colors represent different k 
values with increments 3.6 % of T-Y. (c) Spectra for k = kp normalized as described 
in Ref. [43] and in text. 



found in optical spectroscopy. As with ARPES it is likely that the quasi- Id chains being 
probed in this surface sensitive spectroscopy are protected from surface effects by lying 
well below the cleavage plane. 

Photoemission measurements extended to lower temperatures have continued to 
support LL physics, but have also revealed interesting disagreements with TL-model 
predictions. Lineshapes from angle integrated ARPES data from 300 K down to 30 K 
are very well fitted by TL-lineshapes but only if a is T-dependent [9] . The T-dependence 
of a nicely connects the high T value of ~ 0.9 to the low T-value of ~ 0.6 found in STS. 
As with STS there is no evidence of the resistivity upturn. However, the fact that a is 
T-dependent is clearly outside the one band TL-model. A microscopic theory offered 
in Ref. [9] is that Coulomb interactions involving the two bands of the two chains per 
unit cell (implying four collective modes, two holons, two spinons [31 153|) give rise to a 
T-dependent renormalization of a. 

T-dependent ARPES lineshapes measured over the T range 300 K to 30 K showed 
QC-scaling [43J. However the T-dependence observed was not in agreement with 
expectations from the TL-model. As shown in Fig. l(b-d) of Ref. [IB], although the 
lineshapes sharpened considerably with decreasing T, the sharpening was not as great 
as predicted, even when taking account of the measured T-dependence of a. 

We now summarize aspects of the T-dependent ARPES study that are important 
for the later discussion of Section 5 . Fig. |2]Ja) shows the 300K spectra for wavevector 
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k varying along the T-Y direction of the Brillouin zone. Bands A and B approach Ep 
no closer than 0.12 eV. Bands C and D merge and disperse to cross E-p together. For 
the particular k-path shown, the D band is too weak to observe and so its dispersion is 
sketched as a dashed line based on data [7J El [9] from other k-paths. Fig. [2](b) overplots 
the spectra in the range of the box of Fig. EJ^a) to show the dispersing holon peak and 
spinon edge approaching Ep. The need to measure for multiple T values within the 
lifetime of the sample necessitated noticeably poorer statistics for these spectra than 
for the spectra (cited above pi [3QJ HO]) that were compared in detail to the TL theory 
lineshapes. But apart from the poorer statistics the general features of the spectra are 
the same. 

The T-dependent spectra were tested for QC-behavior as follows. If the general 
scaling form of the spectral function holds, then T~ v A(k,u,T) is independent of T 
if k is chosen for each T so that k/T does not change, i.e. k = (the kp spectra) 
or A; = cT where c is a constant. As described in Ref. [33] the kp spectra were 
normalized to one another by matching the leading edges (0.25 to 0.4 eV) of the B 
band peak, which lies relatively far from Ep and has no apparent T dependence. Fig. 
W^c) shows the normalized spectra vertically offset for clarity. The spectral intensitives 
at Ep were then matched with a multiplicative T-dependent factor whose inverse has 
a power law dependence from which t] can be deduced as shown in Fig. |3]^a). The 
deviations between the data and the power law fit have no systematic pattern and are 
consistent with the effects of the noise in the spectra (see Fig. EJjb)) giving uncertainty 
to the normalization process leading to Fig. EJ^c). When the kp T-scaled spectra are 
plotted vs. (E — Ep)/kp,T the spinon edges align. The scaling behavior of the spectra 
can be visualized from Fig. E^b) which shows for k = the result of a phenomenology 
that is presented in Section 5 and that gives a good description of the data of Ref. |43j . 
We see that the spinon edges scale, but not the holon peaks. Unsealed k-dependent 
spectra were normalized and T-scaled by exactly the same factors as already found for 
the kp spectra and also showed the same general scaling result. In addition to the lack 
of scaling of the holon peaks, one sees from Fig. Et^a) that rj has a value that is essentially 
a rather than the value (a - 1) that would be expected from the TL model [33]. The 
difference between the fitted power law 0.56 and the measured value of a is not deemed 
significant in view of the uncertainty in the normalization process described above. On 
the one hand the observed scaling behavior confirms a generic property of LL physics, 
but on the other hand the deviations from expectations in the TL model raise difficult 
questions for Id theory [TJ. 

This section closes by calling attention to another very important evidence of LL 
behavior that has recently been reported [33]. A pioneering measurement of the ratio of 
the thermal and electrical Hall conductivities in the T range 300 K to 25 K observed a 
spectacular violation of the Wiedemann-Franz law. As explained in Ref. [33] the spin- 
charge separation of LL physics offers a very natural explanation for this failing because 
heat transfer proceeds by both spin and charge degrees of freedom whereas electric 
current involves charge alone. This work is notable for being clearly bulk sensitive and 
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Figure 3. (a) Log-log plot of T-dependence of inverse of scaling factors (circles) needed 
to obtain Figs. 3(b) and (d) in [43], showing power law nearly T a . (b) The T-scaling 
of fcp spectra in the phenomenological description after the intensity is multiplied by 
a factor of l/T 6 . Direct comparison to data is given in Fig.[5]as discussed in text. 

for demonstrating a general transport technique capable of distinguishing between Fermi 
liquid and LL behaviors. 

4. New Photoemission Data Just Above T S c 

Up to now the only single particle spectroscopy data for temperatures not far above 
Tsc — 1-9 K has come from the STS measurement cited above [52] for which the lowest 
T was 5 K. Our previous T-dependent photoemssion studies, angle integrated spectra 
shown in Fig. 2 of Ref. [9], and ARPES spectra shown in Fig. 1 of Ref. [13], extended 
from 300K down to only 30K. In this section we present and discuss two new sets of 
photoemission data taken at temperatures comparable to those of the STS data. These 
new data provide a much higher resolution view of the FS than we have previously 
published in Fig. 1 of Ref. [7] and also allow a more direct comparison to the low T 
STS results. We note that it is very challenging to achieve very low temperatures with 
ARPES because of the complexities of sample holders with multiple mechanical degrees 
of freedom, and the difficulties of providing thermal shielding of the sample while also 
allowing for photon excitation and electron collection. 

The new data were taken with two different ARPES setups where the lowest 
temperatures were T=4 K and T=5 K respectively. One setup is situated in the 
laboratory of the Department of Material Physics, Graduate School of Engineering 
Science of the University of Osaka. Here we obtained angle integrated data with a photon 
energy of 8.4 eV and a resolution of 5meV . The other setup is the MERLIN Beamline 
4.0.3 at the Advanced Light Source (ALS) synchroton. Here we measured with hz/=30 
eV and a resolution of 12 meV. For both experiments, single-crystal Li .gMo6Oi7 samples 
were grown using the temperature gradient flux method [ID] . The instrumentation in 
the Osaka laboratory consists of an MB Scientific MBS T-l microwave excited rare 
gas lamp monochromatized by use of ionic crystals (CaF2 for Kr and sapphire for Xe) 
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|55j . Here Xe was used for the lamp which gives a photon energy of hu=8AeV. Other 
instrumentation consists a low temperature closed cycle He cryostat manipulator and 
a Scienta SES2002 electron energy analyzer driven by an MBS A-l power supply. The 
beamline at the ALS utilizes a low temperature 6-axis sample manipulator cooled with 
an open-cycle He flow cryostat and a Scienta R8000 kinetic energy analyzer. The 
MERLIN beamline has a elliptically polarized undulator which allows one to choose 
arbitrary polarizations of the incident light. Due to the longer elastic escape depth of 
lower kinetic energy electrons the bulk sensitivity of data from the Osaka measurements 
with hz/=8.4 eV is expected to be higher than that for the synchrotron measurements 
with hz/=30 eV. The drawback of that enhanced bulk-sensitivity is a lower cross-section 
of the photoemission process, which requires longer acquisition times compared to the 
synchrotron in order to achieving sufficiently good statistics. In the present case this 
difficulty was exacerbated by the small size of the sample relative to the photon spot size 
of the laboratory setup and the acquisition time difference was ~24 hours compared to 1 
hour. On the other hand, the temperature for the laboratory measurement is essentially 
the same as the lowest T of the STS measurement while its resolution is even higher 
and would be capable of detecting the ~ 8 meV mean field gap of a transition at T m j n . 

We first present the new ARPES results obtained at the ALS. The sample was 
orientated such that the 1-d chains (b-direction) were along the angular axis of the 
detector, which is along T-Y in Fig. H] (c) and was cleaved in situ on the cold cryostat 
in a vacuum better than 8xl0 -11 Torr. The photon polarization was chosen with the 
electric field vector perpendicular to the b-direction. The temperature was maintained 
at T=6 K. A three-dimensional spectroscopic data set I(k x , k y , E) was obtained by 
rotating the sample around the polar-direction. The FS map shown in Fig. |4] (c) was 
produced by integrating the photoemission signal ±6 meV around Ep, essentially the 
width of the energy resolution. The k-resolution is 0.01 A" 1 along the analyzer slit 
(r~Y) and 0.03 A -1 perpendicular to the slit (along T-X), much better than for our 
previously published FS map [TJ. Black means high intensity in the map and we see two 
vertical black lines which represent the FS. In contrast to the theoretical calculations 
by Popovic et al. [6], we see at this resolution no splitting of the FS. Also it appears 
that the FS is essentially straight at this resolution. A Lorentzian peak-fit estimates 
the value of along the FS to vary by not more than 0.006 A -1 , which is less than the 
k-space resolution and less than the variation of the LDA calculation. 

In order to compare with the STS data and for comparison to the Osaka data 
presented below, we angle-integrated the ALS ARPES data. The result is shown in Fig.H] 
(a). One sees the excellent statistics of these data. A line- fit of the spectrum is indicated 
by the solid line through the data points. This fit was made using the k-integrated 
spectral weight of a spin-rotational invariant TL-model [15] with ^=2, consistent with 
previous results described above. The theoretical spectrum was broadened by the 
experimental energy resolution of 12 meV. In panel (b) of Fig. H]one sees the normalized 
residual of the fit, i.e. the difference between the fit and the actual spectrum divided 
by the actual spectrum. The residual is essentially smooth. The a value of ~ 0.65 is 
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Figure 4. (a) Angle integrated photoemission spectra of Lithium purple bronze for 
T=4 K, 6K and 30 K. The spectra at T=4K and 30 K were taken with a resolution of 
5 meV with hv=8A eV. For reference a gold spectrum with the same settings is also 
shown. The angle integrated spectra for T=6 K stems from the same data-set as the 
Fermi surface depicted in (c). This data-set was taken with energy resolution of 12 
meV and hz^=30 eV. All the spectra are generally well fit by the TL model lineshape. 
The gold spectrum fits well with a essentially zero, corresponding mathematically to 
a Fermi edge. Panel (b) shows the normalized residuals for each fit. Panel (c) shows 
the Fermi surface produced by integrating ARPES spectra ± 6 meV around the Fermi 
energy. The green box in the upper left represents the FWHM of the resolution function 
in k-space. The black color means higher intensity. One sees clearly a very straight 
Fermi surface represented by the two vertical black lines. The upper BZ shows the 
Fermi surface as calculated by Popovic[6 . Compared to the theoretical FS, there is no 
sign of a splitting or warping of the experimental Fermi surface. 

comparable to that obtained at 30 K in our previous photoemission experiments [9l [43] 
As found with STS [52] there is no significant difference for T above and below T m i n at 
the resolution of this measurement. 

We discuss the Osaka results next. The samples were oriented by x-ray diffraction 
with the chains (b-direction) along the angular axis of the analyzer. They were cleaved 
in situ on the cold cryostat manipulator in a vacuum better than 2xl0 -10 Torr. As 
the acceptance angle of the spectrometer is ± 7 degree corresponding to ± 0.11 A -1 at 
/iz/=10 eV, and the Brillouin-zone extends along T-Y to 0.57 A -1 in the chain direction, 
corresponding to ~ 37 degree, we measured at angles of 0, 14 and 28 degree with 
the same conditions of time, photon flux and temperature, and afterwards added the 
spectra to obtain angle integration. After finishing the measurements the Fermi-energy 
was determined by evaporating gold onto the sample. The spectra of LiPB for T=4 
K and 30 K are shown in Fig. H] (a) together with the corresponding gold spectrum at 
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T=4 K. One sees the considerably poorer statistics of the LiPB Osaka data relative to 
that of the ALS data or the gold spectrum. Nonetheless the sharp contrast between the 
reduced weight near in the LiPB spectra and the Fermi edge of the gold spectrum 
can clearly be seen. Solid lines through the data indicate line-fits of the spectra using 
the same TL-model as for the ALS k-integrated spectrum. The theoretical spectrum was 
broadened by the experimental energy resolution of 5 meV which was determined from 
the gold spectrum. The a value of the fit is written beside each spectrum together with 
the standard-deviation resulting from the x 2 -fit. For the gold spectrum, we expected a 
Fermi edge and therefore a=0. This small deviation from zero and also the peak above 
Ep in the residual shown in Fig. IHb) might be seen as indicating some systematic error. 

In panel (b) of Fig. H] one sees the normalized residuals of the fits. The numbers 
next to the curves indicate how much the curves were shifted in the vertical direction for 
improved clarity of presentation. Comparison of the residuals again shows clearly the 
poorer statistics of the Osaka LiPB data. Both the T=4 K and the T=30K spectra show 
a deviation from the TL-lineshape at around 5 meV binding energy. This deviation is 
larger for T=4 K than for T=30 K. Since the resolution of 5 meV is better here than 
in the STS experiment it is logically possible that such a feature could have escaped 
detection in the STS spectrum. But a skeptic might point out that because it appears 
also in the 30K spectrum it does not correlate directly with the T-dependent transport 
properties. A skeptic might also claim a hint of this feature even in the gold spectrum 
residual. Thus the question of whether this deviation is real or is due a small systematic 
error, i.e. not perfectly linear behavior of the detector or a contribution from the 
sample holder because the photon spot is large, or is simply the result of insufficient 
acquisition time to obtain better statistics, cannot be absolutely answered. We present 
the spectra as a further confirmation of the overall large energy scale power law behavior 
and regard them as ambiguous on the possibility of structure at the 5 meV energy scale 
of the resolution. We take it as significant that all the new LiPB data sets presented 
here give essentially the same a values ~ 0.65 to 0.7. This is very much consistent with 
the range of values found in STS and perhaps consistent with a slight low T upturn of 
a found previously, as shown in Fig. [3|a) of Ref. [9]. 

5. Phenomenological Description of T-dependent ARPES Data 

The findings of Ref. [13] show that Lio^MogOij is a QC system, but with important 
differences relative to expectations from the one-band TL-model, specifically in the 
exponent of the temperature prefactor and in the lack of the full sharpening predicted 
for decreasing T, which is part of the non-scaling of the holon peaks. 

These differences can be partially described by a phenomenological momentum 
broadening of the TL spectral function, but the required broadening greatly exceeds the 
experimental momentum resolution, as discussed below. The fact that the experimental 
scaling prefactor is T a rather than T*"" 1 ) is a basic motivation to try an integration of 
the theoretical TL spectral function because that will draw a factor of T outside the 



PES and the Unusually Robust Id Physics of L^gMo^On 



13 



integral. Convolving the theoretical spectral function (for which rj = (a - 1)) with a 
momentum window function R(p/p Q ), where po is a T-independent width, and using the 
change of variables p = (recall that v is a constant with units of velocity), gives 



A test (k,u,T) = T c 



R 



(1) 



A constant factor of \jv has been drawn into a redefined A. A test (k, u, T) has the 
T a prefactor observed experimentally but the rest of the expression is no longer a 
universal function of % because T also enters in the form po 



v -f and A 



test 



T a /(*, It can be noted in passing that this exercise also shows the way in 

which a fixed experimental resolution has an increasingly large effect on a QC spectrum 
as temperature decreases. 

Choosing R to be a simple normalized Gaussian with p = 0.065 A -1 and choosing 
the low temperature value of a— 0.6, consistent with the experimental prefactor, yields 
a reasonable description of the experimental data below the binding energy for which 
the intensity from band B becomes important. Fig. [3]^b) shows the general scaling 
behavior of A test , including the experimental energy resolution. After the intensity is 
multiplied by 1/T a and the energy axis is scaled by fceT, the theory curves at kv have 
their edges fall on each other while the peak parts do not, exactly what is observed in 
experiment. The slight deviation from the scaled edge at low temperature above E-p is 
seen in the data of Ref. [33], and as shown there, is due to the experimental resolutions. 
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Figure 5. Fitting the ARPES data of Ref. 43 to A test . (a) Fitting for k =k F at 
various temperatures, (b) and (c) Fitting to both k = k F and k away from k F at 300 
K and 100 K respectively. Parameters are: ^=2 and a= 0.6, v s — 1.9 eVA and po= 
0.065 A -1 , same for all panels. Unit of k is A -1 . 
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Fig. E(a) compares the kp data of Fig. [2^d) for 10 values of T with the 
phenomenological A test and Figs. 5(b,c) make the comparison for four values of k, 
T-scaled away from kp at two temperatures. In these figures, the amplitudes of the 
experimental curves have been adjusted to remove the slight (maximum 15%) prefactor 
differences from T a seen in Fig. [3]^a) and ascribed in the Section 3 discussion of that 
figure to uncertainty in the normalization process. By eye there is generally good 
agreement in these comparisons until reaching the binding energy at which the higher 
lying bands begin to contribute intensity in the experimental spectra. From the spectra 
of Fig. [2](a) and taking account of the width of band B this energy is 0.1 eV to 0.12 eV. 
For k- values away from kp the binding energy of band B increases and so the deviation 
also moves to higher binding energies, nearly 0.2 eV. In spite of this agreement, it is 
important to point out that this phenomenology does not produce the T-dependence 
of a that is measured from T-dependent k-integrated spectra. In fact it works only 
for k-values rather near to kp whereas a much larger range of k contributes to the 
k-integrated spectrum that determines a. Thus the phenomenology is not a substitute 
for a satisfactory microscopic model. 

What is the meaning of this phenomenology? The simplest possibility is of course 
that the momentum is indeed broadened. One could argue that the QC scaling predicted 
in the TL model is being observed except for being modified by the experimental 
resolution. However, the value of po is 5 times larger than the instrumental k resolution, 
computed from the instrumental angle resolution as 0.013 A -1 . Such a large extra 
broadening would have to be ascribed to sample surface quality. For example our 
relatively large spot size might be illuminating regions with slightly different orientations 
on the sample surface. The size of the discrepancy makes this interpretation more 
difficult to defend, but nonetheless the phenomenology shows that the issue of finite 
momentum resolution is very important. ARPES scaling data might evolve with 
progress in sample growth and surface preparation techniques. 

Another possibility is that the phenomenology tries to catch some of the intrinsic 
physics. For example it is clear from the T-dependence of a that the one band TL 
model is missing interactions of importance and we know from band theory and the 
ARPES data that there are two bands crossing Ep, associated with the two chains per 
unit cell. Perhaps the phenomenology can be a guide for future theory. It would be 
very desirable to test the data against the spectral function for a two band TL model if 
it were available. 

6. Robust Id Physics 

The combination of data from ARPES and STS along with the Wiedeman-Franz law 
violation cited above firmly establish generic signatures of LL physics in LiPB for T at 
least down to T min . How does the system achieve the 3d character that must obtain in 
the SC state? At this point it is not possible to make a unified interpretation of all the 
relevant data. One can only speculate in a general way, guided by current Id theory. 
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Recent discussions [251 ESI EE] focus on the idea that the crossover to 3d occurs 
before T$c is reached, either that the resistivity rise is due to CDW formation [24"l 26J, 
which is certainly of a 3d nature, or that a CDW does not occur [251 E3 EE] and that 
crossover is more subtle in the resistivity data [27] and in the mechanism whereby it 
occurs [251 EE]. Here we speculate that Id physics may be more robust, with crossover 
to 3d occurring only with the transition to SC, i.e. that above Tsc there is non Fermi 
liquid behavior rooted somehow in Id with no single particle gap or density wave. In this 
connection one can take note of recent arguments [4j for unconventional SC, possibly 
spin-triplet, based on the finding that the critical field is much larger than the Pauli 
limit. 

Because the Id physics of LiPB appears to lie outside standard Id theory [1] we 
can only try to identify generic possibilities that may be robust for this line of thinking. 
Looking first to single particle spectroscopy, the power laws observed for temperatures 
only slightly above Tsc in the STS and ARPES data with 12 meV resolution are an 
initial temptation to consider this hypothesis. As discussed in Section 4, the 4 K angle 
integrated spectra are ambiguious on the possibility of spectral structure on a lower 
energy scale of the 5 meV resolution, but we note that even this energy scale corresponds 
to a temperature of 58 K, much higher than the T of the measurement or the T min of 
the resisitivity upturn. The energy scale can perhaps be pushed lower by the 6 K 
optical spectroscopy [IT] performed down to 1 meV, corresponding to 11.6 K. This work 
found no single particle gap and no change associated with the resistivity upturn but 
did employ a Drude description of the infrared data. A more recent study [22] , albeit a 
measurement down to T=10 K and with a lower energy limit of 6 meV, also found no gap 
but in addition showed that the infrared behavior is actually non-Drude. Taken together 
these two studies suggest no single particle gap and non-Fermi liquid behavior down to 
an energy scale of 1 meV, but definitive low temperature spectroscopic measurements 
probing to energy scales well below 1 meV are clearly needed. 

The resistivity upturn, with the possibility of a gap less than 1 meV, as obtained in 
Ref. [2E] and shown in Fig. []Ja), is the primary transport evidence suggesting crossover 
due to a small single particle gap. However the highly novel gap function and the T- 
independent constant dc magnetic susceptibility measured down to 2 K are cause to 
consider an alternative description of the upturn. 

Power laws are characteristic features of the QC nature of LL physics. Indeed a 
power law is expected for the resistivity of a LL system. Fig. [T^a) shows that the data 
of Ref. [26J can be described by two power laws, one with positive exponent for the 
metallic part [57] and one with negative exponent for the upturn. Fig. \Bjo) shows the 
T-dependent gap that results from recasting this fit as an Arrhenius law. The gap so 
obtained is very similar to that obtained in Ref. [2E] and the difference is very likely due 
to the difficulty of accurately digitizing the upturn in the resistivity data. In this view 
the unusual gap function is the result of forcing an Arrhenius description onto a power 
law. In particular one notes that at low T the gap must roll over and then actually 
decrease in order to map the very fast rising exponential onto the slower rising power 
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law. 

Such a two-power-law description of LiPB resistivity was first put forth in Ref . [25] , 
where the two power laws were ascribed to two independent LL's, one for each of the two 
bands crossing E-p with each having different values of a and quite different roles in the 
crossover physics. This interpretation neglects the Coulomb interaction that couples the 
chains to give symmetric and antisymmetric holon and spinon modes, as discussed in 
Ref. [3] and Ref. [53] and applied to LiPB in Ref. [9] . An alternate LL interpretation for 
the power law upturn is the old proposal of Anderson localization due to disorder, but 
put in the context of LL physics as a low T crossover of the sign of the resistivity power 
law exponent due to disorder [58J . A great concern of early discussions [121 ELU f!5] 
was that disorder and localization could be incompatible with SC. But, as suggested in 
Ref. [15] , the increased sensitivity to disorder in Id could make it possible for the energy 
scale of the disorder to be much less than Tsc- 

As seems true of all hypotheses concerning the resistivity upturn, this one also 
leaves some aspects of the data unexplained, why the same T min would be found for 
the resistivity along all three axes, why hydrostatic pressure suppresses the upturn (and 
also enhances the SC) [18], and why the upturn can be suppressed |26j by a magnetic 
field. For these aspects of the data one sees the attraction of postulating a gap that 
can be affected by pressure or a magnetic field. Whatever is the proper understanding 
of the upturn, a case can be made that by virtue of a power law dependence on T 
it signifies a continuation of Id behavior rather than a loss. The highly directional 
nature of the magnetic field suppression of the upturn [26] emphasizes the importance 
of the Id character in this temperature range and in this connection one notes also the 
remarkable recent observation [26J that for a non-superconducting crystal a sufficiently 
large magnetic field applied specifically along the Id axis appeared to restore the SC 
with a transition temperature considerably higher than 1.9 K. 

To return to the theme of the opening paragraph of the paper, the greatest issue 
to be confronted in any line of thinking that emphasizes Id physics down to Tsc f° r 
LiPB is the role and the energy scale of t±. The data of Fig. H^c) show that the FS is 
unsplit and straight to a greater extent than predicted in current LDA calculations. If 
one accepts the current lack of evidence for a CDW, and that is the stance of this paper, 
then there is a need to understand why the good nesting of the FS does not result in a 
CDW having all the standard properties that are readily detected for CDWs in related 
materials. Whatever is the mechanism it would apply down to T,sc- Within present 
theory the situation is very difficult but perhaps not quite impossible. We now discuss 
three possibilities. 

In a one band LL model it is well known that the easy route to avoiding crossover 
due to t± is a gap in the spinon mode. The gap energy must be overcome for 
single particle hopping between the chains, which renders t± to be irrelevant in the 
renormalization group sense [3]. But pair hopping that is second order in t± remains 
relevant and can lead to SC [3]. If the energy scale of the spin gap is greater than Tsc 
there can be a transition directly from the spin-gapped LL (a Luther-Emery liquid) to 
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a SC Such a gap in LiPB would seem to be precluded by the T-independent dc 
magnetic susceptibility which has been measured down to 2 K [22J. However, the two 
chain model with four modes, as mentioned in the preceding paragraph, has both a 
symmetric and an antisymmetric spin mode. The symmetric mode that is probed in 
the dc susceptibility can be ungapped for some parameter ranges of this model [53J. 
One could think of the possibility of a gap in the antisymmetric spin mode that has not 
been probed in any transport measurement to date and is small enough to have escaped 
spectroscopic detection but is nonetheless larger than Tsc- Such a gap would also be a 
barrier against single particle hopping between the chains and so might render t± to be 
irrelevant. 

A second possibility is the so-called sliding Luttinger liquid (SLL) [60, EU [62]. This 
is a lattice model for an array of coupled chains and specifically includes a mechanism 
for suppressing a purely electronic CDW along the chains. The physical idea [HU E2] is 
that Coulomb interactions between the chains lead to transverse incommensurate CDW 
fluctuations that modulate the charge density on the chains, and therefore the Fermi 
wavevector and therefore the longitudinal CDW q-vector, which frustrates locking of 
the longitudinal CDW. In this special regime the T=0 ground state retains the LL 
properties of a single chain and is stable against t±, against CDW fluctuations and SC 
fluctuations. A spin gap is not required but of course aids the SLL stability. But in any 
case the parameters of the model must be tuned to bring the system to be near such a 
transverse CDW instability, raising doubt as to whether the model could apply to any 
real material. For other parameter values the stable ground state is a CDW, a novel SC, 
or a Femi liquid. In addition there are difficulties specific to its application to LiPB. The 
model does not address the possibility of a thermal phase transition from the SLL state 
to a SC. Also, as T goes to zero for the SLL state, the resistivity is predicted to become 
infinite transverse to the chains and to become zero parallel to the chains. For the 
transverse direction the resistivity upturn is then consistent with this prediction but for 
the parallel direction one must invoke some other mechanism such as disorder. But one 
must then explain why T m j„ is the same for all directions and one must take account of 
how the transverse interactions may change the effect of disorder [62]. Nonetheless the 
concept that the chain CDW can be suppressed by frustration arising from competition 
between parallel and transverse CDW fluctuations may be more robust than the details 
of the model calculations that have been done and so the SLL remains interesting for 
LiPB. 

In the absence of a spin gap or a scenario like the SLL, the only possibility for the 
line of thinking in this section would seem to be that the energy scale of t± is less than 
Tsc- I n that case there is presently no conclusive theory for or against direct crossover 
of LL behavior to SC. The question entails the interplay of the pair tunneling terms 
needed for SC and the extent to which Coulomb interactions might suppress or enhance 
the pair tunneling. Could the energy scale of t± be this small in LiPB? Current LDA 
calculations P [63] suggest values of t± ~ 30 meV that would preclude this possibility. 
However Id fluctuations on the chains are known to produce a suppression of t± [2ll3~tl63] 



PES and the Unusually Robust Id Physics of L%Q$Mo§On 



18 



to an effective value t±(t±/t) a ^ 1 ~ a \ where t is the interchain hopping. For t w 800 
meV suggested by LDA bands for LiPB and the measured low T value of a ~ 0.6, 
one obtains the very small effective hopping value of 0.22 meV, essentially the same as 
Tsc- Although this estimate is terribly sensitive to the value of a, which is not known 
with such great precision, nonetheless it shows that the effective hopping could be very 
greatly suppressed, consistent with the FS data in Fig. Etc). 

To pursue this idea our current ARPES research aims at quantifying the single 
particle inter-chain hopping and setting a more precise bound on the extent of warping 
and splitting of the FS through yet higher resolution measurements just above Tsc- 
New LDA calculations using the NMTO method [65] and downfolding can be used to 
characterize the t± hopping in greater detail [66]. It is anticipated that the theory 
i_i_'s can then be adjusted in a realistic way to describe the ARPES data and that the 
adjustments required can be compared to expectations from Id theory [67]. Thereby it 
is hoped to obtain a more precise picture of the electronic structure that is giving rise 
to the remarkably robust Id properties and perhaps unconventional SC of LiPB. 
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